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By explicit calculation of the two-loop QCD corrections we show that for singlet axial and vector currents the full 
off-shell (VVA) correlation function in the limit of massless fermions is proportional to the one-loop result, when 
calculated in the MS scheme. By the same finite renormalization which is needed to make the one-loop anomaly 
exact to all orders, we arrive at the conclusion that two-loop corrections are absent altogether, for the complete 
correlator not only its anomalous part. In accordance with the one-loop nature of the (VVA) correlator, one 
possible amplitude, which seems to be missing by accident at the one-loop level, also does not show up at the 
two-loop level. 



1. INTRODUCTION 

Recently, Vainshtein [1] found an important new 
relation between form factors of the (VVA) cor- 
relator matching to all orders in perturbation 
theory, in some kinematic limit, a transversal 
amplitude to the anomalous longitudinal one, 
which is known to be subject to the Adler- 
Bardeen non-renormalization theorem [2]. Later 
Knecht et al. [3] were confirming this kind of 
non-renormalization theorem. These recent in- 
vestigations came up in connection with problems 
in calculating the leading hadronic effects in the 
electroweak two-loop contributions to the muon 
anomalous magnetic moment a M [4-7]. 

The diagrams which yield the leading cor- 
rections are those including a VVA triangular 
fermion-loop (VVA while VVV — ) asso- 
ciated with a Z boson exchange 




and a fermion of flavor / gives a potentially large 
contribution, up to UV singular terms which will 
cancel [81, 
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where a is the fine structure constant, G M the 
Fermi constant, T^f the 3rd component of the 
weak isospin, Q f the charge and N c f the color fac- 
tor, 1 for leptons, 3 for quarks. The mass mf is 
if uif < and m/ if mt > . C / denotes 
constant terms. Since, as granted in the Standard 
Model of elementary particles, anomaly cancella- 
tion by lepton-quark duality J2f N c fQfTsf = 
is at work, only the sums over complete lepton- 
quark families yield meaningful results relevant 
to physics. In any case the quark contributions 
have to be taken into account, and treating them 
as free fermions the leading large log ~ InMz 
drops in sum over each family due to the anomaly 
cancellation condition of the SM. 

However, quarks cannot be treated perturba- 
tive and we expect substantial strong interaction 
effects. A framework to investigate the latter is to 
consider the general structure of the VVA three 
point function 

x(0|T{^(x 1 )V r „(a:2)A p (0)}|Q) (2) 
of the flavor and color diagonal fermion currents 
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where tp is a quark field. 

To leading order the correlator of interest is 
associated with the one-loop triangle diagram 
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plus its crossed (q±, /i <-> q%, v) partner. 

For the static low energy quantity a M = | (g — 
2) p = Fm(0), given by the Pauli form factor at 
zero momentum transfer, the VVA correlator is 
required in the limit 

(4) 
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with 



Indeed, in this kinematic region, the leading 
strong interaction effects may be parametrized 
by two VVA amplitudes, a longitudinal wl{Q 2 ) 
and a transversal wt(Q 2 ) one as functions of 
Q 2 = —q 2 , which contribute as [6,7] 
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where A is a cutoff to be taken to oo at the end. 
Vainshtein [1] has shown that in the chiral limit 
the relation 



wt(Q 2 ) p qcd| 
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which was known to hold at one-loop [9], is 
valid actually to all orders of perturbative QCD. 
Vainshtein's theorem follows from the symmetry 
(p, q <-> fi,q + k) (k — ► 0). Formally, discarding 
regularization problems, the asymptotic symme- 
try derives from the fact that 75 may be moved 
from the A p vertex to the vertex by anticom- 
muting it an even number of times [1]. Thus for 
the quarks the non-renormalization theorem valid 
beyond pQCD for the anomalous amplitude wl 

.,l-loop/- ri 2\ ^N c T q Q q 
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carries over to the perturbative part of the 
transversal amplitude. Thus in the chiral limit 
the perturbative QPM result for wt is exact. 
This may be somewhat puzzling, since in low 
energy effective QCD, which encodes the non- 
perturbative strong interaction effects, this kind 
of term seems to be absent. The term is recovered 
however by taking into account all relevant terms 
in the operator product expansion [3,7]. 

In Vainshtein's kinematic limit, what mat- 
ters is the derivative with respect to k taken at 
k = 0. In this case actually the vertex problem 
reduces to a propagator type problem. In the 
calculation described below we have extended 
this to a genuine vertex type statement at the 
two-loop level. As the extensions of the Adler- 
Bardeen non-renormalization theorem for the 
anomalous Ward identity (VVdA) turn out to 
play an important role in new phenomenological 
applications, we will study in the following such 
possible generalizations by an explicit calculation 
of the leading QCD corrections to the 77Z trian- 
gle. 



The vector currents are strictly conserved 
0, while the axial vector current satis- 
fies a PCAC relation plus the anomaly d^A^ = 
2im ^ l5 ijj + ^£ llU p a F^{x)F^{x). We will be 
mainly interested in the properties of strongly 
interacting quark flavor currents in perturbative 
QCD. Our notation closely follows [3]. 
The Ward identities restrict the general covari- 
ant decomposition of W^vpiqi, 32) into invariant 
functions to four terms 
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with — 53 = q\ + <72 and transverse tensors given 
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The longitudinal part is entirely fixed by the 
anomaly, 

(a) 



/ 2 2 2\ 

wi ( 9i , 9 2 , 9 3 ) = 



which is exact to all orders of perturbation theory, 
the famous Adler-Bardeen non-renormalization 
theorem. The Vainshtein relation is obtained in 
the limit upon identifying 

w L (Q 2 ) = w L (q 2 ,0,q 2 ), 

w T (Q 2 ) ee 4 + V,0,9 2 ) + 4~V,0,g 2 ), 

with Q 2 = —q 2 . 

2. THE CALCULATION 

Here we report on recent progress we made in ex- 
tending non-renormalization phenomenon at the 
two-loop level [10] . For details and further refer- 
ences we refer to the latter paper in the following. 
We perform the calculation with conventional di- 
mensional regularization in d — 4 — 2e dimensions 
and use a linear covariant gauge with arbitrary 
gauge parameter throughout the calculation. Our 
procedure of treating 75 is similar to the one used 
in [11]. We write down all fermion loops starting 
with the axial-vector vertex, and then perform 
Feynman integrals and Dirac algebra without as- 
suming any property of 75 at all. In this way all 
diagrams will be expressed in terms of traces of 
10 combinations of 7 matrices. The prescription 
is sufficient to enable us to arrive at amplitudes 
which have finite limits as d — > 4 in the corre- 
sponding covariant decomposition. After this the 
usual formulas 

Tr[757a7/87/it7i'] = 4teaPnv , Tr [757079] = 



valid in d = 4 dimensions were used. In our con- 
vention £0123 = +1 and (1 — 75)/2 projects to 
left-handed fermion fields. 

Tensor integrals were expressed in terms of in- 
tegrals with different shifts of the space-time di- 
mension [12]. All scalar integrals could be re- 
duced to 6 master integrals by using the Grobner 
basis technique proposed in [13]. The expres- 
sions for the individual diagrams are sums over 
21 terms which are combinations of the 6 master 
integrals 
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multiplied by ratios of polynomials in q 2 and d. 



Here D 1 = kf, D 2 = fef, D 3 



(h - 9i) 2 , 



D± = (fci+g 2 ) 2 , D 5 = (fci-fc 2 ) 2 , As = (fc 2 -9i) 2 
and D7 = (k 2 + q 2 ) 2 . The integrals © form a 
complete set of master integrals needed for the 
calculation of massless vertex diagrams with pla- 
nar topology. The planar integral with 6 denom- 
inators 



A(9 2 ,9 2 2 ,9 3 2 ) 



d d k 1 d d k 2 
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can be reduced to integrals © using 



q 2 eP 6 ^(l-2e)I 3 (qlqlqi)I^(q3- 
+e{P 5 {ql ql q 2 2 ) + P 5 (q 2 2 , q 2 3 , 9i)) 
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i?i satisfies the system of differential equations: 

{x(l - x)d 2 - y 2 d 2 y + [7 - (a + (3 + l)x]d x 
-2xyd x d y -{a + l)yd y - a(3}R 1 = 0, 
{2/(1 - y)d 2 y - x 2 d 2 x + [7 - (a + 13 + l)y]d v 
-2xyd x d y ~(a + (3 + l)xd x - = 0, 

where (3 — 7 = 2a = 2s and 



x = — , y 
% 



(£ g = 9_ g = d_ 

q 2 ' x dx' v dy' 



The general solution of this system may be writ- 
ten in terms of Appell functions F4 

(-q£) 2e Ri = AF 4 {s,2s,2s,2s;x,y) 

+BF±{1 - s, 1, 2 - 2s, 2s; x, y)x 1 ~ 2e 
+CF 4 (1 - s, 1, 2e, 2 - 2e; x, jy)?/ 1 " 25 
+L>F 4 (2 - 3e, 2 - 2e, 2 - 2e, 2 - 2e; z, y), 

where 
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The F4 functions can be simplified yielding 
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(A-Q 3 ) £ VA-Q 3 

where A = \/A and 

A = 1 + x 2 + y 2 - 2xy - 2x - 2y, 
Qi = y+l-x, Q 2 = x+l-y, Q 3 = x+y-1, 

G{z) = 2 iq {s,l~s,2-2s,z), 
F(z, u) = F 1 (1, 1 - s, 1 - s, 1 + e; «, w) , 
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The Gauss' hypergeometric function has an in- 
tegral representation 



G(z) 



T(2-2e) f 1 



du 



r2(i- £ )7 [u{i - u)(i - zu)Y 



(12) 



which is convenient for performing the s expan- 
sion. For the expansion of F{x, y) the relation 

(l + s)(l-y)(l-x)(x-y)F(x,y) = 
(x-y)[l + (l-x)(l-y)e]<j>(x,y) 
+ x(x-y- x 2 y + x 2 )d x <p(x, y) 
+ y(x-y -y 2 + xy 2 )d y (j)(x,y), (13) 

may be used to express this function in terms of 
another F\ function which is more suitable for e 
expansion 

4>{x, y) = Fi (1, -£, -£, 2 + £; x, y) = 

(l + £) / du[{l - u)(l - xu)(l - yu)] s . (14) 
Jo 

An expression for R2 may be obtained by differen- 
tiating the one given for R\ . The hypergeometric 
representation for P5 is obtained by solving the 
first order difference equation with respect to d. 
Details of these calculations will be given in a 
separate publication. Series expansion in e for 
various master integrals to the order needed in 
our calculations was given in [14, 15]. Recently, 
further terms of the e expansion for these master 
integrals were calculated in [16]. 

The sum of all diagrams turns out to be gauge 
parameter independent. In the Feynman gauge 
at 53 = and for arbitrary d, the results of our 
calculation are in agreement with the ones pre- 
sented in [11] diagram by diagram. 
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By applying the prescription outlined above for 
the evaluation of individual diagrams, and us- 
ing Schouten's identities, we are able to reshuffle 
terms to match the tensor structures introduced 
in (J7J, exhibiting manifestly the vector current 
conservation. 

3. RESULTS AND DISCUSSION 

Including one- and two-loop contributions, we 
may represent the form-factors in the form 

(±)/ 2 2 2\ at (±)/ 2 2 2\ 

w y T '(<?!, t? 2 ,g 3 ) = rif N c w\^(q 1 ,q 2 ,q 3 ) 
+ an f N c C 2 (R) W$ far? , q\ , q\) 

Wh{ql,ql,ql) = n s N c w i,L(ql,ql,ql) 

+ an f N c C 2 (R)w 2 , L (qlqlq 2 3 ) (15) 

where a — a s /(47r) includes the QCD coupling 
a s , n,f is the number of flavors and N c the num- 
ber of colors. The quarks are in the fundamen- 
tal representation R and the corresponding group 
theory factor is given by 

C 2 (R)I = R a R a , C 2 (R) = 4/3 for QCD . (16) 

We have been working in the MS renormaliza- 
tion scheme. The singlet axial current = A p 
is non-trivially renormalized because of the ax- 
ial anomaly. It is known [17] that in addition 
to the standard ultraviolet renormalization con- 
stant Zy^g which reads •Z'jjg = 1 in our case (as 
^ms"— 1 — 0(a 2 )), one has to apply a finite renor- 
malization constant Z$ such that renormalized 
and bare currents are related as: 

(J 5 p ) r = Z 5 Z MS (J 5 p ) . (17) 

The counterterms coming from the wave function 
renormalization of quarks and ultraviolet renor- 
malization of the axial and vector currents cancel. 
The finite renormalization constant is known at 
the three-loop level [18]. For our calculations we 
must take 

Z 5 = l-AC 2 {R)a . (18) 

The result of the two-loop calculation after 
adding all diagrams is surprisingly simple and, 



normalized according to ljl5|l . is given by 

2 / 2 2 2\ o 

Q3 w 2,L(q 1 ,q 2 ,q3) = -8, 

w { 2 ~T(qhqi,qi) = -v%t (ii^^i), (19) 

qi& 2 w { -^(ql,ql,qi) =8{x-y)A 

+ 8(x-y)(6xy + A)$W(x,y) 

- A[18xy + 6x 2 - 6x + (1 + x + y)A)]L x 
+ 4[18xy + 6y 2 - 6y + (1 + x + y)A)]L y 

q 2 A 2 w^(ql qlq 2 ) = +8A 

+ 8[6xy + (x + y) A]$ (1) (a:, y) 

- 4[6x + A](x - y - 1)L X 

+ 4[6y + A\(x-y + l)L y (20) 

with L x = lnx, L y — hiy. Expression for 
may be found in [19]: 

(x, y) = j { 2 (Li a (-px) + Li 2 {-py)) 

+ l n V In l±W + \n(px) \n{py) + ^ ) , (21) 
x 1 + px 3 J 

where 

p(x,y) = 2/(1- x-y + X), 

The comparison with the results of the one-loop 
calculation reveals that 

W^vp(qx, q2) \two-ioop — 

i Qfl) |one— loop 

(22) 

Multiplying the sum of one- and two-loop terms 
by the finite factor Z§ we arrive at 

WW (51,92) = WW(<?l,g2) \one-loop (23) 

This is the non-renormalization theorem for the 
full off shell correlator at two-loops. While our 
calculation confirms the relations derived in [3] 
and the non-renormalization theorem @ found 
in [1,7], these findings are not sufficient to explain 
our result valid for generic momenta. 

Taking into account the rather non-trivial mo- 
mentum dependence of the form-factors it is very 
tempting to suggest that it could hold to all or- 
ders of perturbation theory because of the topo- 
logical nature of the anomaly, for example. 
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We would like to stress that the surprising 
relation could be discovered only by keeping the 
general non-trivial momentum dependence. The 
anomalous three point correlator exhibits an un- 
usually simple structure, while contributions from 
individual diagrams are very unwieldy. One can 
expect similar effects for other anomalous corre- 
lators. Since at the order considered the QCD 
calculation is essentially a QED calculation, it 
is highly non-trivial whether this carries over to 
higher orders. For the (VVdA) anomalous cor- 
relator a large number of two-loop calculations 
have been performed, mainly in QED and we 
refer to the comprehensive review by Adler [20] 
and references therein. 

In electroweak SM calculations one would a 
priori expect that renormalizing parameters and 
fields would be sufficient for renormalizing the 
SM. Our calculation shows that on top of the 
standard renormalization, it is mandatory to 
renormalize the anomalous currents Jp by the 
finite renormalization factor Z$ because the lep- 
ton currents and the quark currents pick different 
-Z-factors and if they are not renormalized away 
the anomaly cancellation and hence renormaliz- 
ability obviously would get spoiled. As pointed 
out by Adler and many others [20] the point is 
that there exists a renormalization scheme for 
which the one-loop anomaly is exact. Only in 
this scheme anomaly cancellation and thus renor- 
malizability will carry over to higher orders in 
the SM. Our result shows that due to the neces- 
sity of renormalizing away possible higher order 
contributions from the anomaly also the non- 
anomalous transversal contributions are affected. 
We have shown that at least at two-loops the 
entire contribution gets renormalized away in the 
zero mass limit. 
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